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Simulations with 2+1 flavors of domain wall fermions provide us with the opportunity to compare 
the lattice data directly to the predictions of continuum chiral perturbation theory, up to correc- 
tions from the residual chiral symmetry breaking, m res , and &{a) lattice artefacts, which are rela- 
tively small for domain wall fermions. We present preliminary results for the pseudoscalar meson 
masses and decay constants from partially quenched simulations and examine the next-to-leading 
order chiral extrapolations at small quark masses. The simulations were carried out on two lattice 
volumes : 16 3 x 32 and 24 3 x 64, with the lattice spacing fixed at about 0.1 fm. The subtleties of 
the chiral fits are discussed. We also explore the roles of m res and 0(a) terms in the NLO chiral 
expansions and their effects on the chiral extrapolations for the pseudoscalar masses and decay 
constants. 
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Numerical simulations in lattice QCD are essential for understanding the non-perturbative 
aspects of the Standard Model. However, our ability to obtain the values for physical observables 
is constrained by several factors. One among those is that the input quark masses in the simulations 
are far heavier than the physical values of the up and down quarks, due to the limitations of present 
computing powers. To deal with this requires extrapolations from the simulation results to the 
physical points. Chiral perturbation theory ^ (ChPT) provides a solid theoretical foundation for 
these extrapolations. 

However, uncertainty remains concerning the range of quark masses ChPT formulas can be 
applied to the numerical data, partly because the quark masses in the simulations were too heavy 
until recently. With the emergence of powerful computers, progress has been made towards the 
answer. Simulations with both Wilson fermions [Q] and Staggered fermions [J3j] have shown consis- 
tency between the numerical data and ChPT. Both fermion formulations have limitations, however. 
Wilson fermions have large chiral symmetry breaking effects at finite lattice spacing, and many 
new parameters need to be introduced to account for these effects. Staggered fermions even break 
flavor symmetry, and they require the application of staggered chiral perturbation theory to guide 
the extrapolations. (For a recent review, see Ref. [Q].) In both cases the extra parameters needed 
complicate the fits and may obscure some systematic effects. 

The domain wall fermion (DWF) and the overlap fermion formulations provide the cleanest 
way to probe the chiral limit for various quantities. They both preserve flavor symmetry at finite 
lattice spacing, and the breaking of chiral symmetry for domain wall fermions is exponentially 
suppressed by L s , the size of the extent in the fifth dimension. The amount of the residual chiral 
symmetry breaking, denoted as m res , can be easily measured in the numerical simulations. In this 
proceedings we will show some progress in the chiral extrapolations for the pseudoscalar meson 
masses and decay constants from a series of 2+1 flavor domain wall fermion simulations. 

1. Chiral Expansions for Domain Wall Fermions 

To see how the chiral extrapolations for domain wall fermions should be taken, we need to 
understand how m res enters the chiral expansions. The axial Ward identity for DWF []|] is 

A^(x)0(y)} = 2m f (Jl(x)0(y)) + 2(Jl q (x)0(y)) + i(8 a O(y)}. (1.1) 

where (x) is the (partially) conserved axial current. The term containing 7| ? involves fields in 
the 5th dimension, and measures the overlap between the quark states with opposite chiralities that 
are bound to the two walls at the boundaries of the 5th dimension. As the lattice spacing a goes to 
0, we expect to recover the continuum theory. Thus for low energy Green's functions, 7| must be 
related to the pseudoscalar density 7| by 

Jl CJ = m ies J a 5 +0(a). (1.2) 

With this identity, the terms associated with the input quark mass m f and the residual mass m Tes are 
indistinguishable. Only the "effective" quark mass mf + m Tes plays a role. We can then write down 
the Symanzik effective action to 6 {a) for domain wall fermions: 

S eff = [d\y(x)(W-m q )y(x)]+m- aL °c dwf y(x)oV v Fp V \i/(x) (1.3) 
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where m q = m/ + m res . This is quite similar to Wilson fermions [Q], but the coefficient of the 
dimension-5 Pauli term is exponentially suppressed by L s for domain wall fermions. 

To construct the chiral expansions as done in [Qj, we need to do the power counting properly. 
By dimensional analysis, a rough estimate for the size of the Pauli term is 

ae~ aLs c dwf ~ m res (aA QC D) 2 ~ 0.1m res , (1.4) 

assuming a is about 0. 1 fm and Aqcd ~ 500 MeV. Thus as long as m res is not extremely large, which 
it is not in our simulations, the Pauli term is a higher order effect compared to the mass contribution. 
To leading order, only those terms of €?(m q ) appear in the chiral Lagrangian. However, to next-to- 
leading order, both terms of ^(m q ) and £?(ae~ aLs ) contribute, and the chiral Lagrangian becomes 

^nlo = yTr[£p + (Lp) f ]+J?° LO (1.5) 

where the superscript a indicates the expression at finite lattice spacing and stands for the form 
in the continuum limit, p is a spurion field coming from the Pauli term in Eq. [L3| and has the same 
chiral transformation as the ordinary mass field. This corrects the continuum chiral formulas for 
the pseudoscalar masses and decay constants in the following way: 

• The pseudoscalar masses pick up an additional constant term which is proportional to ae~ aLs ; 

• The expression for the decay constants does not change. 

In the case of 2+1 flavor sea quarks and degenerate valence quarks, the final expressions, at the 
chiral scale ju, can be written as 

( 48 16 
M 2 PS = p+Xv{l + j2(2L 6 -L 4 )x + -2(2L s -L 5 )xv 

1 { 2Xv ~ Xl ~ Xs , Xv (Xv - Xi)(Xv - Xs) „ , Xv 

t 24i¥ [ lv-x n (Xv~Z^ 7" 

+ (Xv - XjKXv - X.) ( , + log ^ )+ (X, - X,)(Xn - X.) 1 (1 6) 

XV — Xt\ A 4 {Xv — Xr] ) M > 



fps = /{l + ^UX + LsXv) 



f 2 



1 



167T 2 / 2 



to + X,)log^ + ^,og^]}. (1.7) 



2\l 2 2\l 



The subscripts / and s denote the light and strange sea quarks, and V denotes the valence quarks. 
Xi = 2B mu for i = V, I, s, or tj, Xn = \(Xl + 2x,s) and % = \QXl+Xs)- Note that here m i is the 



sum of the input quark mass and m res , as implied by Eq. |L3 

2. Numerical Results 

Preliminary work on the chiral extrapolations for domain wall fermions has been reported in 
[^]. There we suffered from low statistics and a large residual mass, thus no reliable NLO chiral 
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fits could be produced. In the past year, the RBC and UKQCD collaborations have generated large 
ensembles of 2+1 flavor domain wall fermion lattices with 5-dimensional volumes of 16 x 32 x 16 
and 24 3 x 64 x 16. These ensembles were generated using the exact RHMC algorithm [§J, and the 
Iwasaki gauge action with j8 = 2.13 was used. For details of the generation of these ensembles and 
their properties, see Ref. [g]. The lengths of the ensembles are typically 4000 MD time units. The 
inverse lattice spacing is about 1.6 GeV [ |l(i| ] which makes the the spatial sizes of the lattices 
2 fm and 3 fm. With the quark masses in our simulations, the finite size effects in these volumes 
are small compared to the statistical errors. Thus in the following discussions, they are neglected. 

2.1 Measurement Details 

Various source smearing techniques have been used to do measurements on these ensem- 



bles [10]. The results presented here are all obtained from the Coulomb gauge fixed source with 
a spatial extent of 16 3 . Partially quenched measurements with valence masses from 0.01 to 0.04 
were performed for both volumes. An additional lighter mass point of 0.005 was included for the 
24 3 x 64 lattices . The number of measurements for each ensemble is shown in Table [jj 



mi/m s 


am va i 


# measurements [sources] 
16 3 x 32 24 3 x 64 


0.01/0.04 


(0.005) 0.01 0.02 0.03 0.04 


704 [2] 


60 [2] 


0.02/0.04 


(0.005) 0.01 0.02 0.03 0.04 


710 [2] 


53 [2] 


0.03/0.04 


(0.005) 0.01 0.02 0.03 0.04 


704 [2] 


55 [2] 



Table 1: Parameters for the measurements. 
Note that the normalization of the Coulomb gauge fixed wall source is not known analytically. 



In contrast to the conventional method of determining fp$ using a point-like source []12|], we need 
to compute the relative normalization of the source operator to the conserved current. Bearing this 
in mind, we can extract the decay constants by 

(0\A L \PS)(PS\A$\0} 



fps 



07 L 7 W 



M PS 



where the superscript W(L) denotes the wall(local) operators. And 



Ex(M)(x,fW(o)) 



(2.1) 



(2.2) 



This method gives statistically more accurate results than the point-like sources. The statistical 
errors for Mp$ and fp$ measured in these ensembles are typically below 1%. 

2.2 Chiral Fits at NLO 



To perform the chiral extrapolations using Eq. [L6J and |1 .7| , we need to know the value of 
m res in the first place. The ratios of the "mid-point" correlator to the pseudoscalar correlator are 
computed on the lattice for each quark mass. And m res is defined as the value in the chiral limit. 
With L s = 16, m res is about 0.003 for these ensembles [10]. We consider two cases for the chiral 
fits: (I) the G{a) corrections are neglected, i.e., the continuum chiral formulas are applied; (II) the 



full formulas of Eq. 1.6 and Eq. 1.7 are used. They are discussed in turn as follows. 
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2.2.1 Fits to continuum chiral formulas 



Neglecting the constant term in Eq. 1.6, the total number of the unknown parameters is 6. Since 



both Mp S and fp$ depend on Bo and /, they should be included in the chiral fits simultaneously. 
For comparison, we first fit the pseudoscalar masses and decay constants independent of each other 
with the 0.04 valence point excluded from the fits (Figure [j]). The data is described by the NLO 
chiral formulas amazingly well. However, the values of Bo and / obtained from the fits of Mj, s 
and fps differ by a factor of 2. This inconsistency comes from the fact that the lattice data of fp$ 
does not have the curvature predicted by the NLO ChPT The quark masses we are working with 
correspond to pion masses of 390 to 630 MeV, which may be outside the chiral region where NLO 
ChPT is sufficient. It is thus not surprising to see how badly the simultaneous fits to Mj, s and fp$ 
fail, as shown in Figure |2[ 





out simultaneously . The data points at am va i = 0.04 are excluded from the fits. 
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The failure to fit our data to the NLO chiral formulas in this mass range motivated us to do 
partially quenched measurements with a lighter valence mass of 0.005 on the 24 3 x 64 lattices. In 
Figure || simultaneous fits for Mp S and fp$ to the 0.005 to 0.02 mass points are plotted. Although the 
fitted curves do not perfectly represent the data, as indicated by the large # 2 /dof in Table ||, we can 
see the tendency of them being more consistent with the data as the quark masses go lighter. The 
fitted curves miss the heavy mass points badly, which may indicate that the NNLO contributions 
are significant beyond the mass range in the fits. 




Figure 3: Partially quenched NLO chiral fits on the 24 3 x 64 ensembles. The fits to Mp S and fp$ are 
performed simultaneously in the quark mass range of 0.005 to 0.02. 



Fit 


£ 2 /dof 


2Boa 


fa 


2Zy6 L4 


2L 8 -L 5 


u 


U 


P/2B 


I 

II 


9.8(3.6) 
10.5(3.4) 


4.12(10) 
3.58(19) 


0.0555(15) 
0.0550(19) 


0.56(13) 
0.83(16) 


3.14(32) 
5.5(1.1) 


-1.05(34) 
-0.66(49) 


4.95(86) 
5.04(93) 


9.2(3.7) 



Table 2: Values of the low energy constants determined from the partially quenched NLO chiral fits to Mp S 
and fps without and with the 0(a) corrections for the 24 3 x 64 ensembles. The quark masses included in 
the fits are from 0.005 to 0.02. The values of L,'s and p /2Bq have been multiplied by 10 4 . The chiral scale 
H is taken to be 1 GeV. 



2.2.2 6(a) corrections 

Since the quark masses we have measured on the 16 3 x 32 lattices are too heavy to have 
consistent chiral fits, in this section we only consider the data for the 24 3 x 64 ensembles. Including 
0(a) effects in the chiral fits only introduces one more parameter p in Eq. L6. We expect the size 
of p /2Bq to be much smaller than m q as discussed in Section |[ Applying such fits to the data does 
not change the fit quality very much. It only slightly shifts the values of the low energy constants 
determined from the fits. Table || compares the fitted parameters from the fits without and with 
the G(a) corrections, denoted as Fit I and Fit II respectively. As expected, the value of p/2Bo is 
about 0.001, which is 1/3 of m res , or 1/8 of the smallest quark mass (0.005+0.003) in the fits. This is 
consistent with the assumption we make to employ the whole power counting scheme that the 6(a) 
chiral symmetry breaking effect is small compared to the mass term for domain wall fermions. 
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3. Summary and Future Plans 

To summarize, the finite lattice spacing corrections for domain wall fermions are found to be 
small, with our current masses and statistics. It is possible to apply the continuum chiral formulas 
directly in the chiral fits. However, the quark masses in current simulations are still too heavy to 
see a perfect consistency between the predictions of ChPT and the numerical results. To answer the 
question of how light the quark masses should be to have reliable chiral extrapolations, we need to 
perform simulations at even lighter quark masses. Our measurements on the 24 3 x 64 ensembles 
with a valence quark mass of 0.001 are in progress. The RBC and UKQCD collaborations are also 
generating a new 24 3 x 64 ensemble with the light sea quark mass of 0.005. These new data points 
will give us a good opportunity to probe the chiral limit for various quantities. We also plan to 
study the chiral extrapolations at NNLO with the full formulas calculated by Bijnens et al. [|3j|. 
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